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Abstract 



Given a minimum measurable length underlying spacetime, the latter may be effectively re- 
garded as discrete, at scales of order the Planck length. A systematic discretization of continuum 
physics may be effected most efficiently through the umbral deformation. General functionals 
yielding such deformations at the level of solutions are furnished and illustrated, and broad 
features of discrete oscillations and wave propagation are outlined. 
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1 Introduction 



A generic theme of robust arguments [T] has established an expectation of a fundamental minimum 
measurable length in nature, of the order of l Planck = ^/TiGn/c 3 ~ 1.6 ■ 10~ 35 m, the corresponding 
time for which is l Planck /c. The essence of such arguments is the following (in geometrical Planck 
units, in which h, c, and rn planck = ^JKc/Gn are chosen to be 1). In a system or process characterized 
by energy E, no lengths smaller than L can be accessed or measured, where L is the larger of either 
the Schwarzschild horizon radius of the system, ~ E, or, for energies smaller than the Planck mass, 
the Compton wavelength of the aggregate process, ~ 1/E. Since the minimum of max(£, 1/-E") lies 
at the Planck mass, E ~ 1, the smallest measurable distance is widely recognized to amount to 



X i 



Planck ' 



Thus, continuum laws in nature are expected to be deformed, in principle, by modifications 
of 0(l P[anck ). Remarkably, even as something like a fundamental spacetime lattice of spacing a = 
0(1 pianck) i s thus likely to underlie conventional physics, continuous symmetries, such as Galilei or 
Lorentz invariance, can actually survive unbroken such a deformation into discreteness, in a nonlocal, 
umbral realization [2]. 

Umbral calculus, pioneered by Rota and associates in combinatorics contexts [SIS], specifies, in 
principle, how functions of discrete variables representing observables get to "shadow" their contin- 
uum limit (a — > 0) systematically, and effectively preserve Leibniz's chain rule and the Lie Algebras 
of the difference operators which shadow (deform) the standard differential operators of continuum 
physics. (For a review useful to physics, see [5].) 

Nevertheless, while the continuous symmetries and Lie Algebras of umbrally deformed systems 
might remain identical to their continuum limit, the functions of observables themselves are modified, 
in general, often drastically so. Usually, the controlling continuum differential equations of physics 
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are discretized [6]. Then, such discrete difference equations are solved to yield umbral deformations 
of the continuum solutions. However, some complication may be bypassed by umbrally deforming 
the continuum solutions directly. Still, as illustrated below on the simplest case of oscillations and 
wave propagation, the modifications may become technically cumbersome in extracting the actual 
physics involved. 

In this brief note, an explicit deforming functional more tractable than the abstract umbral 
deformation definition is articulated, (I17ll34p . which is suited to transforming continuum solutions. 
It is an integral functional shadowing the Fourier representation. It should be of utility in inferring 
wave disturbance propagation in discrete spacetime lattices, by evaluating the umbral transforms 
of standard continuum physics quantities. Can generic features of discrete spacetime impact wave 
propagation at cosmic distances? We first illustrate the umbral language on harmonic oscillations 
discretely deformed through asymmetric finite differences, A + , for simplicity, so the reader could 
quickly confirm all statements by standard numerical analysis textbook methods (Section 2). We 
then point out how the umbral transform kernel technique proposed may bypass conventional umbral 
calculus complications for more realistic symmetric differences A s , for which basic polynomial sets and 
umbral manipulations are less straightforward (Section 3). We conclude by discussing generic features 
of wave propagation on such lattices at cosmic distances, and the logical possibilities afforded by the 
index of refraction modifications at 0(i_p/anc*)i we finally propose conceivable solitonic applications 
(Section 4). 



2 Overview of the Umbral Correspondence 

For simplicity, consider discrete time, t = 0, a, 2a, na, ... Without loss of generality, broadly fol- 
lowing the summary review of [5], consider first the A + discretization (umbral deformation) of d t , 

Ax(t) = *(* + *)-*(*) , (1) 
a 

and whence of the elementary oscillation equation, x(t) = —x(t), namely 

= x(t + 2a)-2x(t + a) + x(t) = 
cr 

and investigate the periodicity of its solutions. Of course, this can be easily solved directly by the 
textbook Fourier-component Ansatz x(t) oc r*, [7], to yield (1 ± a) t//a . To illustrate the powerful 
systematics of umbral calculus [5J, however, we produce and study the solution in that framework, 
instead. 

The umbral framework considers associative chains of operators, generalizing ordinary continuum 
functions by ultimately acting on a translationally-invariant "Fock vacuum", 1, after manipulations 
to move shift operators to the right and have them absorbed in that vacuum. Using the standard 
Lagrange-Boole shift generator 

T= e a9t , so that Tf(t)-l = f(t + a)T-l = f(t + a), (3) 

the umbral deformation is 

T — 1 

d t ~ A = , (4) 
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t i-> tT~\ (5) 
t n ^ (tT~ l ) n = t{t-a){t-2a)...{t- (n - l)a)T~ n = [t]"T" n , (6) 
so that = 1, and, for n > 0, [0] n = 0. These are called basic polynomials [H [6]. 

A linear combination of monomials (a power series representation of a function) will thus trans- 
form umbrally to the same linear combination of basic polynomials, with the same series coefficients, 
f(t) l— > fi^T' 1 ). All observables in the discretized world are thus such deformation maps of the 
continuum observables, and evaluation of their direct functional form is in order. Below, we will be 
concluding the correspondence by casually eliminating translation operators at the very end, through 
operating on 1, so that /(tT -1 ) • 1. 

The umbral deformation relies on the respective umbral entities obeying operator combinatorics 
identical to their continuum limit (a — > 0), by virtue of obeying the same Heisenberg commutation 
relation, 

[dt,t] = 1= [A^T- 1 ] . (7) 

(Formally, the umbral deformation reflects (unitary) equivalences of the unitary irreducible represen- 
tation of the Heisenberg- Weyl group, provided for by the Stone-von Neumann theorem. Here, these 
equivalences reflect the alternate consistent realizations of all continuum physics structures through 
systematic maps such as the one illustrated here.) 

Thus, e.g., by shift invariance, TAT -1 = A, 

[d t ,t n ] = nt n - 1 ^ [A, [t] n T- n ] = nit] 11 - 1 T^ n , (8) 
so that, ultimately, A[t] n = n[t] n - 1 . 

Likewise [31 E] , 

[t} n T- n [t] m T- m = [t] n * [t} m T~ n - m = [ t ] n+m T- {n+m \ (9) 

and so forth. The right member of the equality is the implicit definition of the product [S] defined 
by dotting on 1, [t] n * [t] m = [t] n+m . 

For commutators of associative operators, the umbrally deformed Leibniz rule holds [H [2] , 

[A, fitT-^gitT' 1 )) = [A, f{tT- l )]g{tT- 1 ) + /(tT" 1 ) [A, gitT' 1 )} , (10) 
ultimately to be dotted onto 1. 

Now note that, in this case, the basic polynomials [t] n are just scaled falling factorials a n (t) n , 

so that \-t\ n = (-) n [t+a(n-l)} n . (Furthermore, [an] n = a n n\ ; for < m < n, [t} m [t-am] n - m = [t] n 
; and for integers < m < n, [am] n = 0. Thus, A m [t] n = [an] m [t] n - m /a rn .) 

The standard umbral exponential is then natural to define as, [3] El E], 

E(\t,a) = e A M ee e xtT " ■ 1 = ^ ±j[t] n = ^(Aa) n ( /0 ) = (1 + \af a , (12) 

n=0 ' n=0 ^ ' 
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the compound interest formula, with the proper continuum limit (a — > 0). Evidently, since A • 1 = 0, 

Ae x[t] = X e x[t] , (13) 
and, as already indicated, one could have solved this equation directly to produce the above E(Xt, a) 

a 

Serviceably, the umbral exponential E happens to be an ordinary exponential, 

\ r -I ln(l + Aa), 

e A W = e-^^ 1 . (14) 

The umbral exponential actually serves as the generating function of the umbral basic polynomi- 
als, 

d 



dX> 



-{l + Xa) t/a 



= [t) n - (15) 

A=0 



Conversely, then, this construction may be reversed, by first solving directly for the umbral eigen- 
function of A, and effectively defining the umbral basic polynomials through the above parametric 
derivatives, in situations where these might be more involved, as in the next section. 

By linearity, the umbral deformation of a power series representation of a function formally 
evaluates to 



fit) 



F(t) = /(tT- 1 ) ■ 1 = / (J^j (1 + Xaf a 



(16) 

A=0 



as a consequence. This may not always be easy to evaluate, but, in fact, the same argument may be 
applied to linear combinations of exponentials, and hence the Fourier representation, instead, 



poo r 00 rlk f 8 \ t ^ a 

F(t) = J drf(r) J ^e- fc (l + tkaf a = (l + a- f(r 



t/a 

(17) 

T=0 



The rightmost equation follows by converting k into d T derivatives and integrating by parts away 
from the resulting delta function. Naturally, it identifies with the above eqn ( fT6i) by the (Fourier) 
identity f(d x )g(x)\x=o = gid x )f(x)\ x= o- It is up to the ingenuity of the particular application of such 
functionals involved to utilize the form whose domain of applicability is best suited for the evaluation 
sought. 

It is also straightforward to check that this umbral transform functional yields 

d t f ^ AF; (18) 

or to evaluate the umbral transform of the Dirac delta function, which amounts to a cardinal sine, 
or sampling, function, 

sin(f(l+t/q)) 

6{t) " Ma + t)) ; (19) 



1 There is an infinity of nonumbral extensions of this solution, however: multiplying this umbral exponential by 
an arbitrary periodic function g(t + a) = g(t) will not be visible to A, and thus will also yield an eigenfunction of A. 
Often, as below, such extra solutions have a nonumbral, vanishing, continuum limit, or an ill-defined one. 
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or to evaluate 

/ = ~ F = e l ' a a^Y{t/a + 1, 1/a), (20) 

an incomplete Gamma function, and so on. (For example, discrete integration of the above cardinal 
sine fucntion leads to the umbral transform of the step function, a hypergeometric function 2 -Pi 
of imaginary argument.) In practical applications, evaluation of umbral transforms of arbitrary 
functions of observables may well be more direct, at the level of solutions, through this deforming 
functional, eqn ( fPTI) . 



For example, one may evaluate in this way the umbral correspondents of trigonometric functions, 

e »W _ e -i[t] e »W _|_ e -»[t] 
S ' m[t] = 2i ' °° s[t] = 2 ' (21) 

so that 

Asinft] = cos[t] , Acos[t] = - sm[t\. (22) 

Thus, the umbral deformation of phase-space rotations, 

x = p, p = — x 1— > Ax(t)=p(t), Ap(t) = — x(t), (23) 
readily yields, by directly deforming continuum solutions, the oscillatory solutions, 

x{t) = x(0) cos[t] + p(0) sin[t], p(t) = p(0) cos[t] - x(0) sm[t]. (24) 

As indicated, the umbral exponential being an ordinary exponential; and, by 

(1 + ia) = VTT^ e iarctan(a) , (25) 
these solutions are seen to actually amount to discrete phase-space spirals, 

x(t) = (l + a 2 )^(x(0)cos(ujt)+p(0)sm(ujt)^, p(t) = (1 + a 2 )^ (p(0) cos(ut) - x(0) sin(a;t)) , (26) 

with a frequency decreased from the continuum value 1 to 

uj = arctan(a)/a < 1 , (27) 
effectively the inverse of the cardinal tangent function. 

That is, for 9 = arctan(a), the spacing of the zeros, period, etc, are scaled up by a factor of 

, nS tan(0) , . 

tanc(0) = — > 1 . (28) 
o 

(For complete periodicity on the time lattice, one further needs return to the origin in an integral 
number of N steps, thus a solution of iV = 2irn/ arctana.) Example: For a = 1, the solutions' 
radius spirals out as 2*/ 2 , while to = 7r/4, and the period is r = 8. 

Note that the umbrally conserved quantity is, 

28 = x(t) * x{t) + p(t) * p(t) = x(0) 2 + p(0) 2 = (1 + a 2 )^ (x{t) 2 + p(t) 2 ^j , (29) 
(AS = 0), with the proper energy as the continuum limit. 
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3 More Symmetric Cases 



Unfortunately, the A + of eqn ([T]) is not time-reversal odd, and thus its square is not time-reversal 
invariant — whence the awkward outspiraling of the solutions of the previous section. (Its time- 
reversal conjugate, A_, would have in-spiraling solutions.) 

Instead, one often chooses half the difference of these two delta operators, i.e., the time-reversal- 
odd umbral deformation, 

T - T- 1 

°< " ^ s ^ 2a ' (30) 
The eigenfunctions of A S E = XE are now two [TO] . 

E± = (Xa ± v / l + (Aa) 2 )* /a ; (31) 



one, E + , going to the exponential in the continuum limit; but the other, E_ ( "nonumbral" ) , simply 
oscillating to zero — an oscillation of infinite frequency. 

Now, since [A s ,2t/(T + T" 1 )] = 1, the basic polynomials, [£]" = (t 2/{T + T" 1 ))" • 1, would be 
harder to evaluate, in principle. Thus, they have been evaluated [B] from A s [£]™ = n[t]™ _1 instead, 

n-1 

[t}: = tH(t + a(n-2k)). (32) 

k=l 



However, according to the general generating function consideration of the previous section, they 
may alternatively be generated more directly from E + , 



[C = -(Aa + yi + (W) 



(33) 



A=0 



(E.g., the reader may easily check that [t\l = (t + a)t(t — a), etc.) Since E_E + = (— )*/ a , E_ generates 
the reflected basis [£]"(— )*/ a+n . In general, finding the eigenfunctions of a given difference operator 
and utilizing them as generating functions of basic polynomial sets may provide shortcuts in effort. 

Rather than using umbral deformations of power series representations for functions, however, 
one may instead infer, mutatis mutandis as in f|T7|) . umbral transforms of Fourier representations, 



r°° r 00 riu , , \ t/ a 

ffl ~ F.{t) = J drf(r)J —e-" k [tka+^l-{kaf) , (34) 

to evaluate umbral deformations for general observables, as well as nonumbral ones relying on E-, 
with a minus sign in the above kernel of the deforming functional. 

Thus, now there are four solutions to the analog of eqn ([2]), 

(A 2 S + l)z(t) = , (35) 
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namely 

x{t) = (±ia ± y/l - a 2 )* 7 ". (36) 

Thus the discrete-time solution set 

x (t) = (-) Nt / a (x(0) cos(ut) + p(0) sin(cut)) , p(t) = (-)™/° ^( ) cos(c^) - x(0) sin(wt)) , (37) 

maps onto itself under time-reversal, for integer parameter N = 0,1. (A s x(0) = (—) N p(0).) This 
eqn ( 1351) only connects even points on the time lattice among themselves, and odd ones among 
themselves. Thus, all even points on the time lattice behave the same for even or odd parameter N. 
However, for the N — 1 solutions, the odd time points hop out of phase by it (reflection with respect 
to the origin in phase space), as they are not dynamically linked to the even points^, a phenomenon 
familiar in lattice gauge theory. 

For N = 0, the frequency is increased over its continuum limit: 

to = arcsin(a)/a > 1. (38) 

For N — 1, the arcsine effectively advances by it and the frequency has an additional component of 
n/a. Thus, these nonumbral solutions collapse to in the continuum limit. 

The conserved energy is more conventional, 

2£ s = x(tf +p(t) 2 . (39) 



This reversal-odd difference operator (1301 is the one to be considered in wave propagation in 
the next section, to avoid presumably unphysical exponential amplitude modulations, growths or 
dwindlings, peculiar to the asymmetric derivative seen in the previous section. 



4 Wave Propagation Outlines 



Given the features of discrete oscillations outlined, simple plane waves in a positive or negative 
direction x would obey an equation of the type 0, [10] , 



(Ai - Af) F s = 0, (40) 

with the symmetric difference operators of the type (130]) on a time lattice with spacing a, and an 
x-lattice of spacing b, respectively, not necessarily such that b = a in some spacetime regions. 

For generic frequency, wavenumber and velocity, the basic right-moving waves e*^* _fcx ) have phase 
velocity 

uj a arcsin(fe) . , 

v ( u > k =T/; — > 41 

k oarcsm a 



2 Actually, as above, if f(t) is a solution of ([55)1 . g{t)f(t) will also be a solution for arbitrary periodic g(t + 2a) = g(t). 
Thus, even though (— )'/ a is one such possible g(t), there are even more solutions with arbitrarily mismatched moduli 
(phase-space radii) and phases between the odd and even sublattices — only their frequencies of rotation need be the 
same. The solution set is four dimensional. 
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that is to say, the effective index of refraction in the discrete medium is (6arcsin(a))/(aarcsin(6)), 
so modified from 1 by 0(l p ). 

Small inhomogeneities of a and b in the fabric of spacetime over large regions could yield interesting 
frequency shifts in the index of refraction, and thus, e.g., whistler waves over cosmic distances. It 
might be worth investigating application of the umbral deformation functional fl34|) on such waves, 
to access long range effects of microscopic qualifications of the type considered. 

A further, more technically challenging application of the umbral tranforms proposed might attain 
significance on nonlinear, solitonic phenomena, such as, e.g., the one-soliton solution of the continuum 
Sine-Gordon equation, 

(d 2 x -d?)f(x,t) = sin(f). (42) 



The corresponding umbral deformation of the equation itself would now also involve a deformed po- 
tential sin(/(x 2 i , t ^-i )) • 1 on the right-hand side, for the A s deformation — and 

sin(/(xT~ 1 , tTF )) ■ 1 for the A + deformation. 



As illustrated here, rather than solving difficult nonlinear difference equations, one may instead 
infer the umbral transform that, e.g., the continuum one-soliton solution maps to, 

F s = J d * dT JP dk arctan ^ me v&^| e - ixp - iTk (ipb+y/l - (pb) 2 y /b (ika+ y/l - (fca) 2 )^. (43) 

For the A + deformation, one would have, instead, 

F + = J°° d X dTd P dk arctan f me &\ e^ tTk (ipb + ^(ika + l)^. (44) 

Closed form evaluation of this integral appears complicated, but it could be plotted numerically and 
the representation could yield qualitative asymptotic insights on the 0(lpi anc k) modifications of such 
umbral solitons. 

Likewise, the analog integrals with the continuum KdV soliton f(x,t) = |sech 2 (^(x — vt)) as 
input, 

F, = jT ^y " sech 2 (£<x - „)) e-^(ip b + V^Wf\ika + VWW)"°, 

(45) 

and, for the A + deformation, 

F+ = pxdrdpdk „ sech2 ^v (x _ vt) ^ e - OT - irt ( !p6+1 )" 1 ( ft . a + 1 )"» i (46) 

could be plotted numerically and compared to the Lax pair integrability machinery of Ref [5] , or the 
results on a variety of discrete KdVs in ref [11]. These questions are left for a forthcoming study. 
The ingenuity of further particular applications of the deforming functionals proposed here is left to 
the reader. 
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DE-AC02-06CH11357. Helpful discussions with D K Sinclair, J Uretsky, and D Fairlie are gratefully ac- 
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